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Abstract 

We construct the most general perturbatively long-range integrable 
spin chain with spins transforming in the fundamental representation 
of gl(N) and open boundary conditions. In addition to the previously 
determined bulk moduli we find a new set of parameters determining 
the reflection phase shift. We also consider finite-size contributions and 
comment on their determination. 



Contents 

1 Introduction [1] 

2 The Spin Chain Model d 

3 Constructing Conserved Charges 4 

4 Asymptotic Bethe Ansatz [5] 

5 Degrees of Freedom 8 

6 Tests of the Bethe Equations 13 

7 Finite Length Contributions 15 

8 Summary 16 



1 Introduction 

During the past few years a long list of evidence has been collected that the spectrum 
of planar anomalous dimensions in M = 4 supersymmetric gauge theory is equivalent to 
the energy spectrum of a certain integrable spin chain |1-3|. The spectrum of this spin 
chain is described efficiently by Bethe equations, see |4| and references therein, at least 
for sufficiently long chains. This spin chain model is special in many respects. Most 
importantly the interactions are genuinely non-local and between multiple sites of the 
chain |2|. 

Our knowledge of such integrable long-range chains is still very limited. The best 
known exception is the Inozemtsev chain which is similar to the Haldane-Shastry chain 
but with elliptic rather than trigonometric dependence on the separation of spins. The 
situation for the Inozemtsev chain is similar to the M = 4 gauge theory chain: The 
asymptotic spectrum (for long chains) is described by Bethe equations while determina- 
tion of equations for the exact (finite-length) spectrum remains a challenge. 

Although the M = 4 gauge theory chain is similar to the Inozemtsev chain |5|, it 
belongs to a more general class of long-range chains with simultaneous interactions be- 
tween more than two spins. In the full Af = 4 gauge theory chain the interactions can 
even change the number of spin sites |6|, but here we shall focus on the su(2) sector 
of the model where the length-fluctuations are frozen out. The setup for such pertur- 
batively long-range 1 chains was outlined in |2|. In the simplest case, it is a chain with 
spins transforming in the fundamental (spin-|) representation of su(2). More degrees 

1 The name long-range spin chain actually refers to the fact that for finite coupling A the interactions 
are indeed of infinite, i.e. long range. However, considering the problem within perturbation theory all 
relevant contributions have finite, i.e. short range. Even though it might therefore be more appropriate 
to speak of short-range chains we stick to the notion of long-range chains for historical reasons. 
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of freedom corresponding to a deformation |7| of the overall magnon scattering phase 
were later discovered in |8|. A full treatment of all deformation parameters in the more 
general setting of gl(N) with fundamental spins and up to interactions of range 6 was 
conducted in |9|. It showed that, independently of the rank N, there are four types of 
deformation parameters: The parameters a^(A) govern deformations of the dispersion 
relation. The parameters /3 rjS (A) correspond to deformations of the overall magnon scat- 
tering phase. The parameters 7 r , s (A) fix the linear combinations of commuting charges in 
terms of some canonical basis. And finally the parameters q jW (A) correspond to similar- 
ity transformations of the integrable system without impact on the spectrum. Here the 
coupling constant A ~ controls the range of the Hamiltonian: A contribution at order 
X e is allowed to have interactions among £ + 2 neighboring sites. Complete integrability 
of this system was initially only a conjecture based on the existence of one conserved 
charge. Later it was shown that the Hamiltonian possesses Yangian symmetry which 
constitutes a formal proof of integrability |10|. 

In this paper we shall consider open spin chains which is a natural generalization of 
the above closed chains. In the original gauge theory setup open boundaries correspond 
most naturally to "quarks" , i.e. to fields transforming in the fundamental (as opposed 
to adjoint) representation of the gauge group |11|. On the string theory side of the 
AdS/CFT correspondence, the quarks are represented by strings ending on D-branes |12|. 
However, even in a gauge theory with adjoint fields only, such as M = 4 SYM, open spin 
chains make an appearance: For example one can turn gauge covariant local operators 
representing open chains into a non-local gauge invariant object by means of a Wilson 
loop |13|. Alternatively an open chain (open string) can end on a determinant-like local 
operator (giant graviton) |14|. 

Here we would like to perform an exhaustive study of open perturbatively long-range 
integrable spin chains with gl(iV) symmetry analogously to the one for closed chains 
in |9|. Our aim is to understand how the long-range interactions can deform the boundary 
conditions as well as what the restrictions for the bulk Hamiltonian are. For practical 
reasons we will not study the most general boundary conditions, but only those which 
preserve manifest gl(N) invariance. Note that in many physical models, such as most 
of the systems discussed above, the bulk symmetry is actually broken by the boundary 
conditions, e.g. gl(N) — > gl(M) x g[(N — M). Therefore our results do not apply directly 
to these models, but we expect that the qualitative picture will be roughly the same as 
for our symmetry-preserving boundary conditions. Indeed previous results on concrete 
open perturbatively long-range chain models |15| confirm this expectation. 

Our procedure is the same as in |9|: We will make a general ansatz for two long- 
range spin chain operators including open boundary terms. By demanding that the two 
commute we obtain a tentatively integrable system. Unfortunately, Yangian symmetry 
is broken by the boundary conditions in nearest-neighbor spin chains even if they are 
integrable. Therefore we cannot use Yangian symmetry to provide a formal proof of inte- 
grability of our boundary terms. We will then perform the coordinate Bethe ansatz and 
obtain asymptotic Bethe equations for long-range open chains. These describe the spec- 
trum of the perturbatively long-range Hamiltonian (if indeed the assumed integrability 
holds). 
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2 The Spin Chain Model 



We consider an open spin chain with spins transforming in the fundamental representa- 
tion of gl(N). A spin chain state of given length L is an element of the tensor product 
space (C N )® L . A basis for such states is given by 

|ai, • • • ,a L ), a k = l,...,N, (2.1) 

where each a k represents one basis vector of C . Each homogeneous Ql(N) invariant local 
operator can be built from permutations of adjacent sites. We write these permutations 
in the form [ CL\ , . . . , CL n ] acting homogeneously on a spin chain of length L by 

L-n+l 

[a 1 ,...,a n ]\b 1 ,...,b L ) = \b u ...,b k ,b >. (2.2) 

k=i 

For example each spin chain state is an eigenstate of the operator [1] with eigenvalue L. 
Note that for a closed chain, the operator [1,2] would have the same property. For an 
open chain its eigenvalue is L — 1 instead as will be explained below. 

We can thus introduce the length operator £ counting the number of spin sites in 
the chain. Furthermore we define the operator B to measure whether the chain has a 
boundary or not, i.e. to vanish on a space of closed spin chain states and to equal 1 on 
a space of open chains. In our notation we have 

C = [1], 

B=[l]-[1,2]. (2.3) 
These operators commute with all local operators. 

Boundary Terms. The new feature of open spin chains as compared to the closed 
chains are boundary terms enlarging the set of building blocks for the Ql(N) invariant 
operators. This stems from the fact that spectator legs may not be dropped by identifying 
certain states since these states act differently on the boundaries as illustrated in Figure 1 

[l,ai + 1, . . . ,a n + 1] ^ [ai, . . . ,a n ] ^ [a 1; . . . , a n , n + 1]. (2.4) 

Curiously these boundary terms can be encoded into homogeneous bulk terms, i.e. struc- 
tures of the form 

[1, ai + 1, . . . , a n + 1] - [ai, . . . , a n ], 

[a 1 ,...,a n ,n + l]-[a 1 ,...,a n ], n > 1, (2.5) 
only act on the boundaries and vanish identically in the bulk. 
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Figure 1: The two interactions contributing to a boundary term act differently only at the 
boundaries. 



Type of structures 


Bulk 


Boundary 


Finite Length 


All 


Number of structures 


R\ - (R- 1)1 + 1 


2(R-1)\ - 1 


E(*-2)! 

Z=3 


R 

i=i 



Table 1: Numbers of the different kinds of building blocks with interaction range R in the 
ansatz for a generic local operator in this problem. 



Finite Length Terms. When spectator legs occur on both sides of the operator we 
find another interesting structure. Operators of the form 

[1, ax + 1, . . . , a n + 1, n + 2] - [1, a\ + 1, . . . , a n + 1] 

-[ax,...,a n ,n + l) + [ax,...,a n ], n > 2, (2.6) 

exclusively act on states of length n and vanish identically on the rest of the space of 
states. We are not going to include these terms into our considerations, i.e. we replace 
terms of the form [1, ax + 1, ■ ■ ■ , a n + l, n + 2] by terms of shorter range according to (2.6) 
at every point of our computation. 

Note however that these structures provide a notation for operators acting on chains 
of one specific length only and are therefore useful for addressing the problem of finite 
length effects in perturbative long-range chains. We analyze this issue in Section 7, 

Parity. We introduce a parity operation acting on fundamental interactions as 

P[ax,..-,a n ] P' 1 = [n - a n + 1, n - a n _i + 1, . . . , n - ax + 1]. (2.7) 

It is useful to classify the interactions according to their parity, e.g. if one is interested 
in a parity conserving model. 



3 Constructing Conserved Charges 

The integrability of a spin chain model is expressed by the existence of an infinite tower 
of conserved charges, all commuting among each other. The idea of a perturbatively 
integrable spin chain is that these commuting charges are expressible as a perturbation 
series in a small parameter A |2 

oo 

Q r = J2X k Q { r k \ r = l,2,.... (3.1) 

k=0 
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Here the Q 2r have maximal interaction range 2r + k, i.e. with each power of A the range 
of the charge increases by one. 2 

It is well known that for open spin chains only half of the integrable charges present 
for closed chains are conserved. It turns out that the odd charges Q,2r+i commute with 
the even charges Q,2s only up to boundary terms. Therefore the tower of commuting 
charges can only be constructed of the even charges Q 2r with an even interaction range 
2r at zeroth order in A. 

Making the most general ansatz for the charges with index 2r, i.e. a linear combination 
of all local operators of length < 2r + k with arbitrary coefficients, we require that 

[Q2r,Q2s] = 0(\ k+1 ), r,s = l,2,... (3.2) 

and solve order by order for the coefficients in the ansatz. In this paper we explicitly 
construct the first two commuting charges <2 2 and <2 4 , where Q 2 is defined to be the 
Hamiltonian of the system. Note that technically this is merely a necessary but not 
a sufficient condition for integrability. For a closed chain, however, the existence of 
one conserved charge Q3 has experimentally turned out to be a sufficient condition for 
integrability |9, 10 1. The reason for this behavior is unclear, but we expect the same to 
be true for open chains. 

The construction is analogous to the one for closed chains in |9|. It is somewhat 
more challenging because the range of Q4 is one step longer than the one of Q3 and 
consequently the ansatz contains many more structures at a given perturbative order. 
We present the result for the Hamiltonian Q2 up to second order in Table 5 at the end 
of the paper. We have also constructed Q 2 and Q 4 up to third order, but the resulting 
expressions are too lengthy (and not enlightening) to be presented here. Our result agrees 
with the one for closed chains |9,10| when projecting out the boundary contributions. 
Note however that some of the closed chain parameters are not present for open chains 
because they are incompatible with the boundary conditions. 

4 Asymptotic Bethe Ansatz 

We now perform the asymptotic Bethe ansatz for the open Ql(N) spin chain. For closed 
chains the general ansatz was presented in |16|. There the idea was to construct eigen- 
states of the Hamiltonian for an infinitely long chain as a superposition of asymptotic 
n-particle states. The zero-particle states yield the vacuum energy while the one-particle 
states determine the dispersion relation E(p). Two-particle states with momenta p and 
q then fix the scattering matrices S(p, q) in flavor space such that one ends up with a 
two-particle eigenstate of the Hamiltonian. Integrability implies that the scattering of 
several particles with many different flavors reduces to the two situations where particles 
of either different or the same flavor scatter with each other. Considering a perturbatively 
long-range Hamiltonian 7i(X) the particles can see each other already at finite distances. 
This leads to a scattering matrix which depends on the distance of the magnons in 
contrast to a sharp change of phase for the zero-order Hamiltonian. 

2 We assume the length of the spin chain to exceed the length of the considered operators. Hence, 
our analysis will not necessarily apply to short chains giving rise to spanning interactions, see Section 7 
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e i{ Pk -p k )L = e »0 L (p fc ) S(p k ,qi) S- 1 (p k ,q i ) S{p k ,qj) S~ 1 {p k ,q j ) e-^h) 




Figure 2: Boundary Bethe equations: A particle of momentum p k scatters with an excitation 
of momentum qj and qt. The particle is reflected at the left boundary where its momentum 
changes to p k and it picks up a boundary scattering phase. Having scattered with the particles 
i and j again, this time with a different momentum pk, the particle acquires another phase 
factor at the right boundary and the momentum changes back p k — ► pj~. After this period of 
2L sites, associated with a phase factor of e l ^ Pk ~ Pk ^ L , the magnon has returned to its original 
position, moving in its original direction. 

Imposing closed periodicity conditions, i.e. the closed Bethe equations, on the in- 
finitely long chain then requires that shifting a particle by L sites yields a phase factor 
e ipL which is to be equal to the product of scattering matrices corresponding to the 
interactions with the other excitations on the chain. 

For the open spin chain these periodicity conditions change: In order for the particle 
to arrive at the same position on the chain, moving in the same direction, it has to be 
shifted by 2L sites. On this way the magnon is reflected at the two boundaries and 
passes all other excitations on the chain twice. Each boundary gives rise to a boundary 
scattering phase such that the difference of these phases is part of the Bethe equations. 
Furthermore the momentum p of the ingoing particle changes to p after reflection. 3 The 
two momenta are the two solutions of the equation E(p) = E(p) related by the reflection 
map. For periodicity we therefore have to require that the phase factor e t<yP ~ p ^ L equals 
the product of scattering matrices S(p,qj) and S(p,qj) with the boundary phase factor 
e 2i0(p) ag illustrated in Figure 2. 

Again, integrability implies that the multi-particle problem reduces to two-particle 
scattering and single particle interactions with the boundaries. Therefore it suffices 
to consider a one-particle state, with a single flavor due to gl(N) invariance, in order 
to determine the boundary scattering phase and the reflection map p(p). The bulk 
scattering matrices S(p, q) can be computed in the same fashion as for the closed spin 
chain. 

To determine the boundary quantities we require the state 

L 

|$ L/R ) = ^$ L/R (x)|x>, (4.1) 

x=l 

3 Note that we have p ^ — p since the considered Hamiltonian does in general not preserve parity. 
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with 

$ L (x) = e ipx + e 2^L( P ) e ^-i) + aL ( x ) ? 

$ R (x) = e v(*-£) + e 2*(ri e *-i-i) + aR ( x ) ; (4.2) 

to be an eigenstate of the Hamiltonian close to either the left or the right boundary. 
That is to say that we project the eigenvalue equation 

ft|$ L/R ) = £(p)|$ L/R ) (4.3) 

on all states within the interaction range of the boundary at current perturbative order 
in A and solve the resulting system of equations. This defines the boundary scattering 
phases (fi L/R (p) as well as the reflection map p{p) as a perturbation series in A. Fur- 
thermore we have to solve for the local boundary parameters a L (x) and a R (x) which, 
however, represent UV physics and will not be of further interest in this paper 

a L (x) = 0(X X ), a R (x) = 0(X L ~ X+1 ). (4.4) 

After determining all physically relevant scattering factors for the problem, we can 
impose the periodicity conditions for the considered spin chain type. We first present 
the Bethe equations for the open qI(N) spin chain and afterwards discuss the different 
kinds of free parameters of the system. Each state is described by a set of Bethe roots 
■u^fc. The label £ — 1, . . . , N — 1 indicates the flavor of the Bethe root whereas the label 
k = 1, . . . ,Ke indexes the set of Bethe roots of flavor £. The main Bethe equation at 
level £ = 1 reads 

- \ L K\ / . \ K 2 i 

X l,k\ TT / u l,k — Uij + I , x \ T-r U\ t k ~ U 2 ,j ~ 3 



Hi) n(^^e XP(2 ^, Mlj)) )n 



x exp(2 i 0( M ,, l) ) ( ) II ( T-T-ll »,,))) II f ' U " " 2 " + 1 

V i,k / 1=1 \ 1 ' k 1 'j ' i=i 



3=1 

L 



3=1 v ' ,J ' 3=1 



Ul,k - U 2 ,j 



(4.5) 



Here we have 



x e,k = x ( u ?,k - f), (4.6) 

and the bar represents the reflection map. For the levels £ — 2, . . . , N — 2 the auxiliary 
Bethe equations take the standard form for gl(N) 

K-l-x i K e K t+1 j 

^ _ tt u e,k - u e-i,j - 2 tt u i,k - ue,j + \ tt u i,k — u e+i,j ~ g 
jJl ue, k - ut-u + | Ulk - Ue j _ i 11 Ulk _ Ue+1J + i 

Kt-x - i K e _ . K e+1 _ j 

x TT U i,k - Uj-lj + g TT ^£,fc ~ ^.j - Z JT ^,fc ~ Uj+lj + g ^ ^ 
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and for the final level t = N — 1 



Kn-2 i Kn-i 
"1 T UN-l,k — u N-2,j — 2 TT U N-l,k — u N-l,j + I 

J=l u N-l,h - UN -2,j + § jJl UN-l,k - UN-lj ~ i 

- . i Kn-i - 

x TT u N-l,k - u N-2,j + 2 TT u N-l,k ~ u N-l,j ~ 1 ^ 

7=1 UN-l,k - u N-2,j - \ jJl UN-l,k - u N-l,j + i 

Note that the parameters u^k do not depend on the coupling A. Therefore the Bethe 
equations are deformed only by means of the rapidity map x, the dressing phase 6 and 
the reflection phase 0. These deformation functions are in turn characterized by the free 
parameters c^(A), /3 r)S (\) and #2s+i(A), respectively which are described below. 



5 Degrees of Freedom 

To decode the role of the different types of free coefficients in this spin chain model it is 
helpful to understand how their numbers increase with the order of A. Constructing the 
conserved charges at order X k we solved the equation 

fy^, Qt l) ] = (5.1) 
;=o 

for the undetermined parameters, i.e. the parameters in and Q^ k \ However, to 
determine the overall number of the characteristic coefficients and their parity, it suffices 
to count the constraints imposed by the homogeneous equation 

[5H w ,Qf ) ] + [H {0 \5Qf ) ] = 0, (5.2) 

since they determine the dimension of the space of solutions. To identify the parity of 
the parameters we simply project out the parity even or odd part of the solutions to 
(5.2) and count the remaining linearly independent structures. These coefficients can be 
classified into certain categories which are described in the following. The numbers of 
free parameters in each category are summarized in Table 2, 



Rapidity Map. The rapidity map x(u) shall be defined implicitly by its inverse 

1=0 

Here the parameter functions a^(A) start at order 0(X i ) 

oo 

a e (\) = ^2\ k af\ (5.4) 

k=e 
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Note that the leading order parameters with even indices otffi are not free but fixed by 
the parameters of the previous orders in A 

af = 0, 

4 2) = -aN 1 '. (5.5) 

We believe this pattern to hold for the higher orders as well even though we do not have 
a proposal for how these parameters are fixed precisely in terms of the lower orders. Note 
that the coefficients a 2 £ have odd parity whereas the parity of the a 2 £+i is even. Hence, 
if the system is to conserve parity, then one must set oil = for all even i. 
The inverse map from the w-plane to the x-plane has the form 



u u 
x(u) = — I — 
K J 2 2 



u 

1=0 



The parameters ai(X) are fixed uniquely by the components of «fc(A) in (5.4). Here 
«o(A) starts at order 0(A) and c^>i(A) at 0(X^ l ^ 2 ^ +1 ). The coefficients cti(\) govern the 
propagation of spin flips in the ferromagnetic vacuum. 



Reflection Map. When a particle is reflected at one of the boundaries, its momentum 
changes from p to p(p). For the first perturbative orders we find 

pip) = - p + 8\a { 1] sin 2 V - + 8A 2 sin 2 V - [af - 4a< 1) a< 1) sin 2 V - - 2(a«) 2 ship) + G(A 3 ). 

2 2 2 (5-7) 

Note that the momentum does not simply reverse its sign but receives nontrivial correc- 
tions at higher orders in A. Due to the parity breaking Hamiltonian this is a necessary 
condition for constructing an eigenstate of 0,2- Therefore the reflection map (5.7) is char- 
acterized by the parity odd coefficients «2r- The in- and outgoing momentum represent 
the two solutions of the equation E(p) = E(p). In our notation the reflection map p — ► p 
is simply given by 

u = -u. (5.8) 

Note that the Bethe equations invert under u^k — * ^,fc whereas they are invariant under 
U£j U£j for each j separately. 

Recalling the definition of the elementary magnon charges for the closed Ql(N) spin 
chain |9 

<£ l0SC » = ( , \ f i " , \. v A , (5-9) 

r — 1 \x{u+^) T 1 x(u — \y v 

we define the elementary magnon charges for the open Ql(N) chain by 

q r (u) = \ (q?° sed (u) + (-l) r q?° scd (-u)) . (5.10) 

This definition is motivated by the transformation behavior under the reflection map 

q 2r {u) = q 2r (-u), 
g 2 r+i(«) = -<?2r+i(-«). (5.11) 
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Dressing Phase. The dressing phase, a common phase factor of the bulk and bound- 
ary scattering matrices, is defined by 

oo oo 

9(ui,u 2 ) = ^ ^ / 9 r , s (A)(g r (ui)g s (u 2 ) - q s (ui)q r (u 2 )) ■ (5.12) 

r=2 s=r+l 

Here the free parameters starting at order (9(A S_1 ) are given by 

oo 

A, S (A) = A *#2- ( 5 - 13 ) 

fc=s-l 

The dressing phase occurs in the open Bethe equations only in form of the difference 

0(ui,k, uij) - 9(u!, k , uxj). (5.14) 

Note that this combination inverts its sign under — > «i,fc whereas it must remain 
invariant for Uij — > Wij. The products of two odd elementary charges (?2r+i?2s+i are n °t 
invariant under U\j — *> U\j and are therefore not allowed 

02r+l,2s+l — 0. (5.15) 

In fact one can convince oneself that these terms are incompatible with the boundary 
Yang-Baxter equation. Furthermore, products of even elementary charges q2r<l2s drop 
out in the combination (5.14), i.e. the parameters /?2r,2s with even index pairs do not 
appear in the Bethe equations. This fact is related to the appearance of bilocal similarity 
transformations which are described below. 

Reflection Phase. The reflection phase as it appears in the above Bethe equations is 
defined by 

oo 

4>{u) = 22 <Wi(A)g 2 s+i W - 0(u, u). (5.16) 

8=1 

The related intrinsic parameters S 2s +iW starting at 0(X S ) are then given by 

oo 

W(A) = (5.17) 

k=s 

These coefficients govern the scattering of particles at the boundaries. Only the difference 
20 of the right and left boundary phase R and L appears as a physical parameter in 
the Bethe equations. Note that the dressing part of the reflection phase can be regarded 
as a j = k contribution to the second line of the i = 1 Bethe equations. Due to the 
transformation rules for the elementary magnon charges we find 

(j)(u) = -<f>{u). (5.18) 
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Eigenvalues of the Charges. If we define the bare charges by 

Q2 S = ^2q2s(ui,k), (5.19) 
the eigenvalues of the spin chain charges are determined by the formula 

oo 

Qir = 72r,- 2 (A)5 + 72ri0 (A)L + ^72r, 2s (A)g 2s - (5.20) 

8=1 

The functions 72r,2s(A) are given by 

oo 

!2r,2sW= E Xk ^r,2s- (5-21) 

k=2 max(s— r,0) 

They correspond to the linear coefficients in the mixing of the bare charges Qis form- 
ing the spin chain charges Qir- Since all bare charges commute with each other, the 
coefficients 72r,2s(A) are not shared between the spin chain charges. All other types of 
coefficients are shared between the charges. 

Note that in contrast to the closed spin chain parity conservation does not impose 
further restrictions on the parameters 72 r ,2s since only the even elementary magnon 
charges qi s contribute to the spin chain charges. The second bare charge Q2 is printed 
up to second order at the end of this paper. 

Local Similarity Transformations. The coefficients and (k,i do not appear in 
the Bethe equations and therefore they do not influence the spectrum. They correspond 
to perturbative bulk or boundary similarity transformations of all operators 

00 

Q r = TQ r T-\ 1=1 + ^^, (5.22) 

k=l 

where is an arbitrary bulk or boundary interaction of range k+1 parametrized by 
or respectively. Contributions to T which are linear combinations of the commuting 
charges do not alter the charges. Thus for counting purposes one has to remove these 
trivial similarity transformations. 

As indicated above, only the even charges Q,2r are conserved quantities for the open 
spin chain. The charges with odd indices Q,2 T +i can be constructed by requiring that 
they commute with the Hamiltonian up to boundary terms. Hence, structures corre- 
sponding to these odd charges do correspond to extra similarity transformations which 
have to be treated as trivial similarity transformations only in the bulk and there have 
to be removed for the counting. For the boundary instead, the odd charges do provide 
nontrivial similarity transformations and have to be added for the counting. 

As an example we can commute the essential structures corresponding to the zero- 
order charges Q2 and Q3 with each other 

Q 2 ~ [2,1], 

Q 3 ~ [3,1,2] -[2,3,1], (5.23) 
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to get 



[Q 2 ,&]~[1,3,2]-[2,1,3]. 



(5.24) 



The resulting boundary structure is not of interaction range 4 as one might expect but 
of range 3. As a consequence Q3 already appears as an extra similarity transformation 
at order A. Similarly the commutator of Q 2 and Q 5 already appears at order A 2 instead 
of A 4 since it is of interaction range 4 instead of the expected range 6. 

Bilocal Similarity Transformations. In order to match up the total numbers of 
parameters it is important to consider a curious class of bilocal similarity transforma- 
tions which exists only for open chains. Bilocal operators can be built from two local 
interactions acting at different positions of the spin chain 

[ai, . . . ,a n |6i, . . .,b m ]\ci, ...,c L ) = 

L—n—m+l L—m+1 



In a natural notation we denote by [Q2r|Q2s] the bilocal composition of the two local 
charges Q 2r and Q 2s where Q 2r acts towards the left of Q 2s . If we commute this particular 
bilocal operator with another local even charge Q 2t we get a local operator. This is due to 
the fact that contributions to the commutator vanish as long as both parts of the bilocal 
operator are well separated. Only if the distance between both parts in the sum above 
becomes less than the interaction range of the local charge, the resulting interactions 
contribute to the commutator. Therefore these special bilocal similarity transformations 
give rise to additional local terms in the commuting charges. The dressing parameters 
@2r,2s with even index pairs do not appear in the Bethe equations and hence do not 
influence the spectrum. 



In the preceding sections we have constructed a tentatively integrable open spin chain 
Hamiltonian and asymptotic Bethe equations to diagonalize it. However, we have no 
rigorous proof for the integrability of the system and thus we cannot be sure that the 
Bethe equations are indeed correct. 

In order to test the Bethe equations (4.5) we have computed some explicit solutions 
for a small number of excitations on short open chains up to second order in A. Plugging 
these Bethe roots into (5.20) we have obtained the energy of the corresponding states. 
On the other hand we have diagonalized the combinatorial Hamiltonian in Table 5 on 
a suitable basis of states and found the spectrum to be in perfect agreement with the 
results from the Bethe equations. 

To be more precise, let us consider a state of length L = 4 with K\ = 2 and K 2 = 1 
Bethe root excitations as an example. It corresponds to the following Young tableau 




/ , | c lj •■) c ki c fc+ai, ••, Cfc+a n , Cfc+n+l, ••, Q, , .., C/ + fe m , Q+, 

l=k+n 



-in- 



.+!,.., c L ). (5.25) 



k=l 



6 Tests of the Bethe Equations 



(6.1) 
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Table 3: Decomposition of tensor products of k fundamental representations of gl(N) into 
irreducibles. 



from the tensor product decomposition of C 4 , cf. Table 3. The level two Bethe root W2,i 
can be eliminated by solving the level two Bethe equation (4.7). We can then solve the 
Bethe equations (4.5) for the Bethe roots ui t i and to find three admissible solutions. 
This is done order by order in A. Substituting the Bethe roots into the energy formula 
(5.20) finally yields the energy of the corresponding state. 

On the other hand we can apply the Hamiltonian Q2 to a basis of states characterized 
by the magnon numbers and flavors we considered when solving the Bethe equations. For 
our example one can take this basis to be given by the states resulting from permutations 
of the entries in |1, 1,2,3). Note that this yields twelve eigenvalues of Q2 as opposed to 
three values we want to compare to. Nine of them, however, represent energies of the 
following Young tableaux in the decomposition of C 4 , cf. Table 3, 

(6.2) 



5 



mm. 



We have compared the energies calculated in these two different ways for all states 
corresponding to the following Young tableaux in the tensor product decomposition 
illustrated in Table 3: 



4cpa 



i, 6cm, 5cj 



rm, 



n. 



(6.3) 



Our tests confirm that the Bethe equations describe the spectrum of Q2 correctly giving 
us confidence that the derived spin chain Hamiltonian is indeed integrable. 
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7 Finite Length Contributions 



Let us discuss briefly the issue of spanning terms (2.6). These span the whole length of 
a finite open chain in analogy to the wrapping terms for closed chains, cf. |17, 18 1. They 
can be represented in the same basis of local operators and thus a complete enumeration 
of such terms is easily feasible by not identifying any of the interaction terms (2.2). We 
have repeated the analysis of Section 3 including the spanning terms. We have found no 
restriction on the bulk or boundary terms, but there are additional spanning terms. We 
will not specify these terms explicitly, but merely present the number of such terms in 
Table 12 

Let us compare the numbers to the set of irreducible representations of Ql(N) for a 
chain of length L in Table 3 (we assume that N > L). We write the decomposition 
formally as 

N L = n Ljfc R iifc , (7.1) 

' 'k 

where riL,k are the multiplicities of the irreducible representations Rl^. Now it turns 
out that the number of irreducible components of the tensor product 

V n Lh = 1,2,4,10,26,76,... for L — 1,2,3,4,5,6,... (7.2) 

^— *k ' 

appears to match precisely with the number of independent coefficients affecting the 
eigenvalues of Q2,4 in Table 4. This means that starting from a certain perturbative 
order the eigenvalue spectrum is freely adjustable. The spanning order for Q 2r at length 
L is given by A L+2 ~ 2r . In other words, our version of perturbative integrability completely 
loses its constraining power at finite length. Equivalent claims have been made in Section 
6.5.4. of |17| for wrapping interactions in closed chains. 

Before we comment on the implications of this result, let us discuss the remaining 
finite length parameters labeled by S. These appear first at order \ L where the spectrum 
of all charges Q 2r becomes fully adjustable. In that case, also similarity transforma- 
tions of the space of chains of length L are permissible. The number of such similarity 
transformations is given by the formula 

Y j n L , k {n L , k -l) = L\-Y j n Ljk = 0,0,2,14,94,... for L = 1,2,3,4,5,... . (7.3) 

'k *— »« 

It is based on the fact that similarity transformations can only act on identical irreducible 
representations. For n identical irreps, there are n 2 similarity transformations, but n of 
them act trivially on a given set of commuting matrices Q 2r - Furthermore the number 
of similarity transformations is equal to the number of spanning interactions of length L 
which equals L\. 

Exceptionally we find a non-zero entry in Table 4 for similarity transformations at 
L = 5 and A 3 . It is explained by a degeneracy of two eigenvalues of Q 2 at leading order 
for identical irreps at L = 5 (the one with multiplicity 6 in Table 3). In addition to 
adjusting these two eigenvalues one can also adjust the eigenvectors with 2(2 — 1) = 2 
similarity transformations. 

How can we understand the complete arbitrariness of finite length contributions? At 
spanning order for open chains (and analogously at wrapping order for closed chains) 
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Table 4: Number of spanning terms at given order and length. Terms listed under Q.2,4 affect the 
eigenvalues of Q2 4j respectively. Terms listed under S correspond to similarity transformations. 

the interactions act on the whole chain. One thus loses the desired locality property of 
Q r and all linear operators become admissible contributions. We may thus consider the 
charges Q r as free general matrices. In the generic case, two n x n matrices commute 
if their eigenvectors coincide. Therefore one should expect 2n degrees of freedom for 
the eigenvalues of the two matrices and n(n — 1) degrees of freedom for the common n 
eigenvectors (up to rescaling). 

What implications does this observation have on integrability of (perturbative) long- 
range chains? There are two possibilities which crucially depend on the precise definition 
of integrability (which is currently unclear): In the one case, integrability permits all the 
finite length terms that we have found. This implies that the operator Q2 is almost 
completely arbitrary at sufficiently high orders of A or at finite A. Clearly there cannot 
be a method to determine the spectrum which is remotely as efficient as the Bethe ansatz. 
Therefore this definition of integrability cannot be useful in practice. 

There may be various stronger definitions of integrability which would constrain some 
or all of the spanning (wrapping) terms. In addition they might even constrain some 
of the bulk (or boundary) terms. The last option would be particularly desirable for 
AdS/CFT integrability because it would provide us with further constraints on the fi- 
nite size spectrum. Consider for example the Haldane-Shastry and Inozemtsev integrable 
long-range chains |19|. Their finite size spectrum is uniquely defined by the Hamiltonian. 
More conveniently it can be described exactly by equations similar to Bethe equations 
(for Haldane-Shastry the equations are purely combinatorial while for Inozemtsev the 
equations are currently known only in special cases). It appears unlikely that the equa- 
tions can be deformed in such a way that only the spectrum at spanning or wrapping 
order is deformed. Hopefully a suitable integrability constraint can be found to constrain 
the spectrum of M = 4 SYM at finite size. 

8 Summary 

In this paper we have generalized the analysis of long-range integrable 0((iV) spin chains 
started in |9| to the case of open boundary conditions. For these open chains only the 
even integrable charges Q2 r are conserved and we have explicitly constructed the first 
two of them up to order A 3 . We have performed the nested asymptotic Bethe ansatz for 
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the open spin chain and found several new features characterizing its moduli space: 



• a new type of deformation parameter f^s+i starting at 0(X S ) which represents a 
phase for reflections of particles at the boundaries of the chain, 

• a reflection map p(p) corresponding to a nontrivial change of particle momentum 
at the boundaries due to a parity breaking Hamiltonian, 

• boundary and bilocal similarity transformations having their origin in the open 
boundary conditions. 

On a set of states we have demonstrated that the Bethe equations yield the correct energy 
spectrum which indicated that the proposed Hamiltonian is indeed integrable. In addi- 
tion we have briefly discussed the role of spanning terms in the context of open chains 
as the counterpart of wrapping interactions for closed spin chains. We found that these 
finite size contributions do not impose any restrictions on the bulk or boundary interac- 
tions but merely provide additional degrees of freedom not influencing the integrability 
of the system. 
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<2 2 (A) = ([1,2] -[2,1]) 

+ ia„(A)([2,3,l]-[3,l,2]) 

+ ttl (A) ([1] - 4[1, 2] + 2[2, 1] + [1, 3, 2] + [2, 1, 3] - [3, 2, 1]) 
+ §<f 3 (A) (2[1] - 2[1, 2] - 2[2, 1] + [1, 3, 2] + [2, 1, 3]) 
+ C 1 ; 1 (A)(-[1,3,2] + [2,1,3]) 

+ 2a (A) 2 ([1] - 2[1, 2] + [1, 3, 2] + [2, 1, 3] - [3, 2, 1] 

+ [2, 3, 4, 1] - [2, 4, 1, 3] - [3, 1, 4, 2] + [4, 1, 2, 3]) 
+ «!(A) 2 (-12[1] + 32[1, 2] - 7[2, 1] - 11[1, 3, 2] - 11 [2, 1, 3] + 8[3, 2, 1] 

+ [1,4, 3, 2] -[2, 3, 4,1] + [2, 4, 1,3] + [3, 1,4, 2] 

+ [3, 2, 1,4] -[4, 1,2, 3] -[4, 2, 3,1]) 
+ |a (A)ai(A) (-18[2, 3, 1] + 18[3, 1, 2] + 3[2, 4, 3, 1] + 3[3, 2, 4, 1] 

-3[4, 1,3,2] -3[4,2, 1,3]) 
+ |5 3 (A) 2 (-2[2, 1] - 2[3, 2, 1] + 3[1, 3, 2] + 3[2, 1, 3] - 2[1, 2, 4, 3] 

-2[2, 1,3, 4] + [1,4, 3, 2] + [3, 2, 1,4]) 
+ 15 5 (X) (-4[1] + 4[1, 2] + 2[2, 1] - 2[3, 2, 1] + [1, 3, 2] + [2, 1, 3] 

- 2[1, 2, 4, 3] - 2[2, 1, 3, 4] + [1, 4, 3, 2] + [3, 2, 1, 4]) 

+ a 1 (X)S 3 (X) (-10[1] + 10[1, 2] + 10[2, 1] - 4[1, 3, 2] - 4[2, 1, 3] - 2[3, 2, 1] 

- [1, 2, 4, 3] + [1, 4, 3, 2] - [2, 1, 3, 4] + [3, 2, 1, 4]) 

+ |/? 23 (A) (-4[1] + 8[2, 1] - 2[2, 3, 1] - 2[3, 1, 2] - 2[2, 1, 4, 3] - 2[2, 3, 4, 1] 

+ 2[2, 4, 1, 3] + 2[3, 1, 4, 2] - 2[3, 4, 1, 2] - 2[4, 1, 2, 3] 

+ [2, 4, 3, 1] + [3, 2, 4, 1] + [4, 1, 3, 2] + [4, 2, 1, 3]) 
+ ie+ 1 (A)([2,4,l,3]-[3,l,4,2]) 

+ *e+ 2 (A) ([2, 4, 3, 1] - [3, 2, 4, 1] - [4, 1, 3, 2] + [4, 2, 1, 3]) 

+ !Ci;i(A) 2 (-2[2, 1] + 3[1, 3, 2] + 3[2, 1, 3] - 2[3, 2, 1] - 2[1, 2, 4, 3] - 2[2, 1,3,4] 

+ [1,4,3,2] + [3,2, 1,4]) 
+ ^(A) (2[1, 2, 4, 3] - 2[2, 1, 3, 4] - [1, 4, 3, 2] + [3, 2, 1, 4]) 
+ < 2 - 2 (A) (2[2, 3, 1] - 2[3, 1, 2] - [1, 3, 4, 2] + [1, 4, 2, 3] - [2, 3, 1, 4] + [3, 1, 2, 4]) 
+ < 2 +(A) ([1, 3, 4, 2] - [1, 4, 2, 3] - [2, 3, 1, 4] + [3, 1, 2, 4]) 
+ «i(A)Ci; 1 (A) (-2[1, 2, 4, 3] + 2[2, 1,3,4]) 

+ 0(X 3 ) 

Table 5: Normalized Hamiltonian printed up to second order. 
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